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Abstract. We study the Dirac equation for spinor wavefunctions minimally coupled 
to an external field, from the perspective of an algebraic system of linear equations 
for the vector potential. By analogy with the method in electromagnetism, which 
has been well-studied, and leads to classical solutions of the Maxwell-Dirac equations, 
O | we set up the formalism for non-abelian gauge symmetry, with the SU(2) group and 

the case of four-spinor doublets. An extended isospin-charge conjugation operator is 
defined, enabling the hermiticity constraint on the gauge potential to be imposed in 
a covariant fashion, and rendering the algebraic system tractable. The outcome is 
an invertible linear equation for the non-abelian vector potential in terms of bispinor 
\Q | current densities. We show that, via application of suitable extended Fierz identities, 

the solution of this system for the non-abelian vector potential is a rational expression 
involving only Pauli scalar and Pauli triplet, Lorentz scalar, vector and axial vector 
current densities, albeit in the non-closed form of a Neumann series. 
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1. Introduction 

One of the central equations of electrodynamics is the Dirac equation, the spinor 
solutions of which describe the states of a free relativistic spin-| particle, such as an 
electron. The electromagnetic vector potential can be introduced by imposing local 
U(l) gauge covariance, and since elements of this group commute, the electromagnetic 
interaction is an abelian gauge field. In order to describe the way in which an electron 
interacts with its own self-consistent electromagnetic field, one approach to solutions 
of the coupled Maxwell-Dirac equations is to "invert" the U(l) gauge covariant Dirac 
equation to solve for in terms of if), then to substitute this into the Maxwell equations, 
which have the Dirac current as the source term pQ. The full Maxwell-Dirac system is 
very complicated, and no exact, closed form solutions have yet been discovered, although 
global solutions have been shown to exist [2J. Solutions for certain simplified cases 
using the algebraic inversion technique have been derived by Radford and Booth, with 
restrictions such as a static Dirac field with further assumptions of spherical (Coulomb) 
[3] and cylindrical (charged wire) symmetry p[J. It is interesting to note further that 
a spherically symmetric external Coulomb field requires the existence of a magnetic 
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monopole term in the vector potential [3J, [5]. Another approach for finding solutions to 
the Maxwell-Dirac system has been pursued by Legg ;(»] . where the algebraic inversion 
was handled without fixing the gauge. The actions of various simply transitive subgroups 
of the Poincare group were then applied in order to reduce the system of equations to a 
set of algebraic relations, involving the group invariant components of a general gauge 
invariant vector field. A global, closed-form solution was discovered for one of these 
subgroups, which exhibited unbounded laminar stream flow of the relativistic electron 
four-current along the z-axis, with mutually perpendicular electric and magnetic fields, 

!■ 

With regards to the inversion of the Dirac equation itself, in an early study [7] 
Eliezer pointed out that, when writing the Dirac equation in the form 

MA^ = R p (1) 

where M contains only Dirac spinor components, because det M — 0, the four equations 
are not linearly independent, and so M is not invertible. However, a more recent study 
by Booth, Legg and Jarvis [8] demonstrated that the system of linear equations is indeed 
invertible if strictly real solutions of the vector potential are required. In the current 
study, we have assumed this condition holds for the 577(2) non-abelian case. 

A major publication by Takabayasi [9] explains that by writing the Dirac equation 
in terms of Lorentz scalars, pseudoscalars, vectors and axial vectors, the model becomes 
analogous to a relativistic hydrodynamical one, and the current solutions describe fields 
of relativistic charged flow. Non-abelian fluid flow was investigated in [10], where an 
application to a quark-gluon plasma, involving the SU(3) gauge group, was discussed. 

In the current study, we are attempting to derive an inverted form of the Dirac 
equation, for interactions involving the SU(2) gauge group. Section 2 contains a brief 
summary of the inversion for the abelian case, outlining the method used in [8], where 
the regular and charge conjugate versions of the Dirac equations are multiplied by 
V'7 / i and ij^ c 'j ll respectively, allowing one to work in terms of bispinor current densities, 
such as ip'j^ip- The rank-2 skew tensors are then eliminated via subtraction and 
the utilization of charge conjugate bispinor identities for the semi-classical case of 
commuting wavefunctions. The well-known inverted form is then easily achieved by 
inverting the remaining Lorentz scalar. An alternative inversion, obtained as a result of 
multiplication by 4 ! l5lfi is a l so given, as well as several consistency conditions. 

Section 3 sets up the inversion of the non-abelian Dirac equation for the SU(2) 
gauge field W afJ ,, and follows the same path of logic as the abelian inversion, where we 
primarily focus on forming bispinors with an extra doublet degree of freedom, such as 
\E'r a 7^\E', by multiplying by ^7^7^. An extended form of the charge conjugate equation 
is required in order make the form of the gauge potential covariant, dubbed the "isospin- 
charge conjugation". This method turns out to yield a much more difficult inversion 
problem compared with the abelian case, involving the expression of an inverse matrix 
of form (J — iV) -1 by a Neumann series. Conditions of convergence are stated, but not 
verified for this case, and a closed form of the inverse matrix has not been obtained 



Algebraic inversion of the Dirac equation for the vector potential in the non-abelian case3 



as of yet. We end the section with a brief comment on alternative bispinor-forming 
multiplication options, which yield couplings of various non-abelian current densities to 
Wan- A full list of coupling equations is deffered to appendix B. 

In section 4 the Fierz expansion of the term is given with the motivation 
of developing Fierz identities, describing relationships between various non-abelian 
current densities. In particular, general Fierz identities for non-abelian Lorentz vector 
current and axial current density products, J ifl Kj U , are derived in appendix C using this 
expansion. Using the abelian case [TT] as a guide, various antisymmetric combinations of 
these identities are then constructed to produce expressions for the rank-2 skew tensor 
current density and its dual *Sj flv , solely in terms of Lorentz scalars and vector 
combinations. Pauli singlet (i = 0) and vector triplet (i = 1,2,3) are treated case-wise, 
due to their different algebraic behaviours. The rank-2 skew tensor current densities 
may then be eliminated from the expression for the inverse matrix (/ — N)^ 1 entirely. 

Conclusions and prospects for future study are given in section 5, followed by a 
robust appendix. The first two appendix sections provide lists of commonly used Pauli 
and Dirac matrix identities, as well as a list of field-current density coupling equations, 
such as those discussed in section 2, but for a greater variety of multiplication options 
and algebraic combinations. The last two sections of the appendix include more detailed 
derivations of the Fierz identities discussed in section 4. 

2. The abelian inversion case 

The Dirac equation for a spin-| fermion of charge q moving under the influence of an 
electromagnetic field is 

(i$ -q4- m)ip = (2) 

which is covariant under a U(l) (abelian) gauge transformation. Here we invoke the 
Feynman slash notation <fi = 7 v cv Rearranging into a more convenient form 

-fipAu = <\> (3) 

where <p = q^ 1 ^ — m)i[), we can obtain an expression in terms of Dirac bispinors by 
(J3]) by -?/>7 M to form bispinors and applying the identity 7^7^ = 8^ — ia^, such that 

W>Ap - fya^Av = (4) 

The form of the charge conjugate Dirac equation is the same, but with the sign of the 
charge in C reversed. Note that we define the charge conjugate spinor as usual as [12] 

^ = CV T = i 7 2 7 °^ T . (5) 

With this in mind, we can define a relationship between bispinors involving charge 
conjugate spinors and regular bispinors 

^T^ c = -tfC-^Cip (6) 
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for any T in the Dirac algebra. In order to perform the inversion, we apply the charge 
conjugation identities 

C-'jJC = -j, (7a) 
C- 1 ^ = -cv (7b) 
C- 1 ( 1 , lu ) T C = lul , (7c) 

to the appropriate bispinor current densities if)Tif), to obtain the relationships between 
the bispinor densities appearing in (jlj) and their charge conjugates, which written 
explicitly are 

^ c = (8 a) 

Note that the allocation of signs would be different if we were to treat the spinors if) 
as field operators if> instead of semi-classical wavefunctions, however the quantum field 
aspects of the Dirac equation are beyond the scope of this paper. Now, subtracting the 
charge conjugate of (BJ from (TJJ itself, applying the identities (|8aj)-(l8dl) and rearranging, 
we obtain the well-known inverted form of the £7(1) covariant Dirac equation 

" 2q W U 

where = if)^^, the probability current density. There is an alternative inverted form, 
which involves multiplying fl3]) by ^757/^ in which case the mass term vanishes 

A = i ^757^ + ^757^ / 1q n 

where 75 = i7o7i7273- By considering spinors that have zero real inner product with 
^ffjip, we can obtain the following consistency conditions [5] 

d,f = (11) 
= -2im^7 5 ^ (12) 

r l5 $ip = = 0, (13) 

where = ipjsjf/if!, the axial (pseudo) current density and ip c = ip T C. These 
consistency conditions constitute the conservation of current and partial conservation 
of axial current, as well as a complex consistency condition originally derived by Eliezer 

®- 

3. Non-abelian SU(2) case 

The 57/(2) gauge covariant Dirac equation for a doublet spinor \l> is 

[i$ - (g/2)r ■ \fr - m]* = (14) 
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where r a (a = 1,2,3) are the non-commutative generators of infinitesimal rotations in 
doublet space and W ail are the Yang-Mills fields, the £77(2) gauge fields analogous to 
Ap. Explicitly, r a are the Pauli matrices, which obey the commutation relations 

[r a /2,r b /2} = ie ab c r c /2 (15) 

where e abc is the rank-3 Levi-Civita tensor, and is antisymmetric under exchange of any 
two indices. Note that throughout this paper, we use the Einstein summation convention 
for the Pauli indices, where a repeated index implies summation, with raising/lowering 
used to highlight the fact. Now, to derive the charge conjugate of this equation, we 
must follow the same process involved in charge conjugating the abelian Dirac equation. 
The goal in the abelian case was to flip the sign of the charge relative to all the other 
terms, and was achieved by complex conjugating the entire equation, then multiplying 
it by an invertible matrix U, such that U'j^U^ 1 = —7^ and Uip* = ip c . We follow the 
same process in the SU (2) case, but not with the presupposition that the sign of g will 
necessarily be flipped. Complex conjugating and rearranging (JUJ), 

{[i<9„ + (g/2)r T ■ W V YT + ™W = 0- (16) 

We have used the fact that r is Hermitian (r* = rj). Multiply f fT6|) by I ® U, with I 
being the 2x2 identity in doublet space, indicating that U commutes with r and acts 
only on the Dirac spinor degree of freedom 

= {[id v + (g/2)r T ■ Wvp'f*!/- 1 + m}im* 
=4- = {[id l/ + (g/2)r T -W U ]Y -m}V c . (17) 

To make the form of the gauge potential covariant, another step is required which is not 
present in the abelian case, to convert the r T back to r. This is done by multiplying 
( ITT)) by e = ir 2 , such that we make use of the Pauli identity r a = —erje^ 1 , 

= {[id„ + (g/2)eT T e- 1 ■ W v Yl v - m}eV c 

0= (g/2)T -ty -m]¥ c (18) 

where we have defined e\l/ c = \I/ 1C as the isospin-charge conjugate (henceforth, IC) spinor. 
Note that the sign of the coupling constant g has reverted back, so that ( TT8)) is of exactly 
the same form as ( TT4j) . Mimicking the abelian case ()3j), we rearrange ( I14j) and ( ITS)) into 
the more convenient forms 

Yt ■ w v y = r b Y^W bu = $ (19) 
Yt ■ w u w c = r b Y^ ic w bu = $ ic (20) 

where $ = 2g~ 1 (i$ — m)^, $ 1C = 2g~ 1 {i$ — m)^ lc and we have replaced the triplet vector 
dot-product notion with a sum over Pauli components b. Multiplying equation (TT9)) by 
^Talfi, applying the Dirac identity that produced (TTJ), as well as the Pauli identity 

r a r b W bll = (6 a b + ie a bc T c )W bfl = W a>M + ie a bc r c W btM (21) 

we obtain the expression 

W aM - Wa^W av + ie a bc ^r c W b , + e a bc ^r c a^W bu = ¥t o7m $ (22) 
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with the form of the IC equation being exactly the same. This is similar to (jl]), but 
with "extra" terms contracted with the Levi-Civita symbol on the left-hand side. We 
require a non-abelian analogue of the bispinor relationship (J6]), which is 

W c {n <g> T)^ ic = -^(e-V^e) g> (C^I^C)* (23) 

with T being an element of the Dirac algebra as before, and r, an element of the Pauli 
algebra, where i = 0, 1,2,3 and i = corresponds to the 2x2 identity. Explicit sign 
relations for particular values of (r^ Cg> T) are not given here, although they can easily be 
calculated following the same method as in ((6]). Following the same method as with the 
abelian case, we subtract the IC of ( 122|) from ff22l . apply the appropriate current sign 
relationships from ( 1231) . then rearrange to get 



^gW at , + e a bc ^r c a^gW bu = i(*r«7 M ^ - #Ka7^) - 2m^r a7 ^. (24) 

We will now define a suite of non-abelian currents, that will be used throughout the rest 
of this paper: 

J l = ^ Ti m (25 a) 

J* = *Ti7 M * (256) 

<V = *T f <v* (25 c) 

'Si^ = *7i75<v¥ (25 d) 

^ = *ri7B7 M * (25 e) 

AT* = *T i75 *- (25/1 
There are 64 current densities altogether, excluding *Si / ^ from the count, since 

*<V = (i/2)e^ p Sr. (26) 

However, a consequence of the results presented in the next section is that the number 
of linearly independent current densities is lower than 64, as the Si^ v terms in any 
expression can be eliminated entirely. Equation ( l24"j) can now be rewritten in the more 
compact form 

(J <W - e a cb S c /) 9W bu = i(*r a7 ^# - *K»7m*) " 2mJ ^- (27) 



It is apparent that there is an additional matrix term on the left-hand side of (1241) that is 
not present in the abelian case, which prevents us from immediately finding an inverted 
form for W atl . Now, let us consider, as with the abelian case, an alternative formulation 
of the above expression, by pre- multiplying ( TTTJ]) and (J2"0j) by ^ r r a 757 At . Applying the 
sign relationships via (123]) and subtracting the IC equation from the non-IC equation 
gives 

(K S^S a b - e a cb *S c , u ) gW bv = i(*r a75 7 M ^ + *K a757 ^), (28) 

in which the mass term vanishes, as in the abelian case. We could go a step further and 
add the two equations ( 1271) and ( 1281) . then divide by the scalar terms to give 

^a Cb {^S c ^ U + *S C/J U ) 



S U S 6 

Jo + K 



bu 
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The matrix on the left-hand side is of the form (J — iV), which is invertible by way of 
the Neumann series 

00 

{I - N)- 1 = ^N n = I + N + N 2 + ... (30) 

n=0 

for iV e C nx ", which has condition of convergence p(N) < 1, where p(N) is the 
spectral radius of N [13J . In the following section, we will show that the terms in 
the expansion of (J — A r )~ 1 can be converted from contractions involving spin-current 
tensors to contractions involving J, K Lorentz scalar and vector currents exclusively, by 
deriving appropriate Fierz identities. Firstly, a brief comment on some of the various 
other objects which the Dirac and IC Dirac equations, (fl9j) and (120]) . can be multiplied by 
to form bispinor-field coupling expressions. We could multiply either of these equations 
on the left by objects of the form ^rjT, where V is an irreducible element of the Dirac 
basis (including the dual of the rank-2 tensor, 7507^), then either add or subtract the 
resulting equations. For example, consider the case where r = J, i = 0, then we multiply 
(fl9|) and (|20|) by ^ and \1/ 1C respectively. Subtracting the IC equation from the non-IC 
equation gives 

(g/2) J bu W bu = (i/2)(*#tt - tf^tf ) - mJ , (31) 
and adding them gives 

TjT^tf = _*j7$#. (32) 

We could choose to eliminate the bispinor with the left-acting derivative operator by 
substituting the second equation into the first, resulting in the expression 

(g/2) J hv W hv = - mJ . (33) 

This equation describes the coupling of the Lorentz vector current J ail with the vector 
potential field via contraction of both the Pauli and Lorentz indices, resulting in a sum 
of Lorentz and Pauli scalar terms. Unlike the invertible cases discussed above, this 
equation has no free indices and therefore can not be inverted via multiplication by an 
appropriate matrix. Since these types of equations may provide valuable information in 
future studies, all multiplication options are listed in appendix B. 



4. Non-abelian Fierz identities 

Consider the 8x8 matrix formed by the product In the SU(2) doublet degree of 
freedom, this is a 2 x 2 matrix, and in the Dirac spinor degree of freedom, it is 4 x 4 
matrix. In the pure Dirac case, the product of two Dirac spinors can be expanded via 
a Fierz expansion 

16 

H> = J>* r « = (V4)(W + (1/4)(^)7 M + (l/SX^vVK" 

R=l 

- (1/4)(^ 757 ^)757 M + (1/4) (^75^)75- (34) 
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There are 16 terms in the sum, and the coefficients an are the Dirac bispinors ipT^ip, 
multiplied by a numerical constant, with Yr being the i?th element of the sixteen- 
component basis of the Dirac algebra. It is interesting to note that extensions of Fierz 
expansions such as this have been described in arbitrary higher dimensions in [13] . Now, 
in the pure 2x2 case, the Fierz expansion for the matrix formed by the product of two 
doublet spinors v is 

3 

vv \ = Y^c % n = (l/2)(^)r + {l/2)(v\ a v)r a (35) 

i=0 

with a = 1, 2, 3, to make four terms in the sum in total. As discussed before, r a are 
the Pauli matrices, and r is the 2x2 identity. Coefficients q are pure SU(2) isospin 
bispinors v^TiV. Now, the basis of the Fierz expansion is the tensor product of the 
Dirac and Pauli bases 

3 16 
t=0 R=l 

+{l/W)S iflv {a^ <g> r l ) - (l/8)i^( 75 7 M ® r*) + (l/8)^( 75 <g> r*) (36) 

where the coefficients are the SU(2) bispinors as previously defined, and are derived 
by pre-multiplying ( 1361) by an element of (r# <g> Tj), utilizing trace identities, and then 
solving for the leftover am. Henceforth, we will exclude the tensor product symbol 
explicitly, however its presence is implied in any product of Dirac and Pauli matrices. 

The Fierz expansion can used to expand products of non-abelian currents in terms 
of other currents in the Dirac-Pauli algebra, for example 

J»K h v = ^7^^757^ 

= -(1/8) J^-f-fTaT^ + (1/8) J^757 M 7 CT 7^arV^ 
-(l/16)^ 757 ^VrarV^ 

+ (l/8)K i(T ^YrYr a T l r b ^ + (l/S^WV^S* (37) 

which, after a lengthy expansion, converting all of the Dirac and Pauli matrix products 
to sums of irreducible terms gives 

J/K b v = (l/A)[U a *Sr + iJ&V - iK a Sr - iK b S/ v + J»K b v 

+J a "K b » + J b »K a » + J b v K a » - JarKfrT - h*K°*r 

+5 ab (U *S ^ - U C *S C ^ - iK Sr + iK c S c » u + WRo" 

+Jo u K > M - J/K cv - J C V K C » - J aK a V^ + JcaK ca n^)} 

+ (l/4)e ab c [-U Q K cV ^ + \K Q J c rT + J c .Jox^ x + K ca K 0X e^ x 

+(l/2)i(-S fl <J *Sf u - SoVS/" + SSSSo™' + SSSSo"")]. (38) 

In the set of Fierz identities for JK vector products, we call this the a — b case. There 
are three other cases of vector current products: J^Ko", J^K a u and Jq^Kq 1 , which 
we call the a — 0, — a and — cases respectively. Due to their long-winded nature, 
we leave the derivation and listing of these Fierz identities to appendix C. Now, if we 
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wish to express the inverted form of the vector potential equation (I2"9"j) independently of 
the rank-2 spin current (skew) tensor S^ u and its dual, we must obtain an identity to 
describe Si^ solely in terms of J, K scalar and vector current densities. Since is 
antisymmetric under exchange of its two Lorentz indices, we should form antisymmetric 
terms from the two- vector product J a fJ, Kb u , then Fierz expand using ( 1361) and solve for 
So"", S a »\ *S ^ or V. 

Let us define the suite of 16 abelian currents as follows 

ct = ^/> (39 a) 

3n = ^1^> (396) 

■V = </"V/'0 ( 39 c ) 

*V = intPiw^ (39cf) 

K = ipl5lv4> (39 e) 

u = ijyy s1 j, (39/) 

where the dual of the rank-2 skew tensor current may be calculated as in ( 126]) . If we 
consider the Fierz identity for the abelian version of the rank-2 spin current tensor [TT] , 

[ae^" K - hj(5/5 v « - 6^6 u ")} 3p k K 



a 2 — lo 2 



(40) 



we can see that, in the non-abelian case, we should consider Fierz expansions of 
antisymmetric current combinations with Lorentz structure of the form J P K U — J U K P 
and e pupK J p K K . Due to the presence of the extra internal Pauli index, we need to 
take consideration of how this will vary the form of Si' 11 ' compared with s pu . As 
discussed in appendix D, the correct approach is to treat the derivation of S pu and S a f ,v 
separately. For the i = case, we calculate the Fierz identities for Jq p Kq u — J a u K ati , 
and J a p K au - J U K P , then add to form 

JfK* - = 2i(J *V - K S ^). (41) 

The antisymmetric part is calculated by adding e pupK Jq p Kq k and e pvpK J ap K a K to form 

^PKj. pK i K = 2 {J Q Sr - K *S n- (42) 
Taking the combination 
J e pupK J ip K\ - iK Q (J^K iv - J»K ip ) 

= 2{jisr - wsr + wsr - K 2 s Q n, m 

the middle two terms on the right-hand side cancel, and we can rearrange to obtain the 
expression 

$T = (1/2)(J 2 - K 2 )- l [J e p % K - iK (S/5 K v - 5 P V 5^)} J^K* K . (44) 

The abelian (1401) and non-abelian (jUJ) cases share a close similarity, with the main 
differences being the factor of 1/2, and the sum over the internal Pauli index in the 
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non-abelian case. Similarly, we can switch the terms that Jo an d K Q multiply in (1431) to 
obtain an expression for the dual 

*S(T = (1/2)(J 2 - Kl)- 1 [K^ v pK - iJ ( W - 5 p v 5/)]JfK l \ (45) 

Note that we can check the validity of this dual identity by using the defining identity 
( |26|) on either (1441) or (145]) . Now taking a step further, the scalar currents inside the 
square brackets can be eliminated by taking the sum or difference of (1441) and f j45l ) 

Sr ± *Sr = (1/2)(J =F #0)" W t KS P V - 5 p ^)]J/K iK . (46) 

Now, in order to form an expression for S a ^ u , we again need to consider antisymmetric 
combinations of JK vector current products, but in such a way as to have the Pauli 
vector triplet index i = a = 1,2,3 present in first order only, combined with the Pauli 
scalar singlet index i = 0. That is, we shall be dealing with the rank-1 JK Pauli vector 
combinations, the — a and a — CclSGS, clS opposed to the — rank-0 or a — b rank-2 
Pauli index cases. As discussed in appendix D, the appropriate Lorentz antisymmetric 
combinations are 

= u *sr + iJa*sr - \K sr - iK a sr m 

as well as 

e^ K {J ap K 0K + J 0p K aK ) = W + J a S<T - K *S a ^ - K a *Sr- (48) 
Taking the combination 

J e^ K (J ap K 0K + J 0p K aK ) - iK Q [WK Q v + J^K a v ) - {J a v K Q » + JfK»)\ 

= (J 2 - Kl)Sr + Wa - K Q K a )Sr + (K J a - J Q K a )*S<r (49) 

we can see that after rearranging to solve for Sa^, we need to substitute the identities 
(|44|) and (145!) to eliminate the Pauli singlet [i = 0) skew tensor current densities. The 
final expression is 



Sr = (J 2 - i^rVo^V - i#o(W - SSSSMJSKo" + J p K a K ) 

'2\2 I 



+ TW^%Yl l K ^ U p« + - S p V)} J/K iK , (50) 



which bears less resemblance to the abelain case ( 1401) than the Pauli singlet case ( 1441) 
does. To calculate the dual, we follow exactly the same process, but switch the Jo and 
K Q in (1491) . which after some rearrangement and substitution gives 

*Sr = (Jo - Kl)- 1 ^^^ - Uo(6 p V - 8 P "8 K »)]WK Q K + Jo p K a K ) 
J 2 + K 2 



+ 2(J 2 -K 2 ) 2 [Ka€flUpK + U *W S " ~ 6 pV)W Km 
- J oK L, 2 [Ja^ ' P . + iKa(6 p »6 K v - 8 p »8 K »)}J/K lK . (51) 
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Again, this can be confirmed by using ( I26p . As in the i — case, taking the sum 
or difference of these two identities removes the current terms from inside the square 
brackets. Following some straightforward algebraic manipulation, we obtain a somewhat 
simpler form 

± v = (i/2)( Jo =f i^r W T - W)l 

•[2( J =F K )(J a p K K + J ?K a K ) - (J a =p K a )Ji P K iK ]. (52) 

The relative simplicity of ( 1521) compared with ( 1501) and ( T5~TT) is the reason we chose to 
pursue an invertible Dirac equation of the form (|29|) as opposed to (127)1 or (|28|) alone. 
In particular, the extra simplicity will have a profound impact on the complexity of the 
higher-power terms in the Neumann series form of the inverse matrix. 

5. Conclusions 

In this study, the inversion of the Dirac equation to solve for the vector potential has 
been extended to the non-abelian case, the 577(2) gauge field specifically. An extension 
of the charge conjugation operation was made in order to make the form of the gauge field 
generators covariant, and the algebraic system tractable. In analogy to previous studies 
performed on the abelian case, the non-abelian Dirac system was re-written in terms of 
bispinor current densities, by way of multiplication by the terms ^^7^ and \t r T a 7s7 /i . 
Combining these equations, an invertible form was achieved, provided we made use of a 
Neumann expansion to describe the form of the inverse matrix. In order to eliminate the 
rank-2 skew tensor current densities from the Neumann expansion, we were motivated to 
derive appropriate Fierz identities by considering antisymmetric combinations of J^Kjy 
current density products. Expressions for S 0/JiU , *S 0flu , S afMU and *S afll/ were subsequently 
derived, and the convenient linear combinations So^i*^^ and S aflu ±*S afJil , were formed. 
Some options for further work include deriving consistency conditions, analogous to 
those presented at the end of section 2, for the SU(2) well as describing 

explicitly the conditions for the convergence of the inverse matrix Neumann series. A 
broader study of the non-abelian Fierz identities is also in order, in particular obtaining 
a complete minimal set, from which all other redundant Fierz identities can be derived. 
The scope of the inversion may also be extended to S77(2) x £7(1) and S77(3) gauge 
fields. In the latter case, the lack of an extended analogue of the SU (2) isospin-charge 
conjugate spinors \1/ 1C suggests that the analysis of Fierz identities will have to confront 
a yet more involved set of bispinor current densities, and the Dirac equation, a more 
complicated inversion calculation. 

Appendix A. Pauli and Dirac Identities 

Here follows various Pauli and Dirac matrix identities that are used throughout this 
paper. 
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Appendix A.l. Pauli Identities 



T a n = S ab + ie abd r d (A.l) 

T a T c n — r a 5 bc + r b 5 ac - r c 5 ab - ie abc (A. 2) 
Appendix A. 2. Dirac Identities 

yyy = ^ 7 a + v ux^ _ _ i^ 757ff ( A _ 4 ) 

yy/V = V^ u V at + V ua V^ - - ^ u o- ut - ir] ua a pt 

+ir ) » a a l/e + \tf^a av + i^V 7 + irfW - ie^ ae ^ 5 (A.5) 

rf a ^ = \rf^ v - ir^-f + e^ e<T 757 CT (A.6) 

a^rf = _ [ v ^y + e^ e<T 7 5 7 CT (A.7) 

= vrfrT - \rf v rf e + - rf v a^ - e aeu ^ 5 + ie aeuX ^a p x (A.8) 

- e xpae e x ^ T = rf»rf v rf T - r) p »rf v rf T + r] pu r] aT r] tp - r] pu r] tT r] ap + rfrf^rf" 

-rTrfrT (A.9) 

Appendix B. Current-Field Coupling Identities 

Here we list the system of expressions that result when we multiply the SU(2) gauge 
covariant Dirac equation and its isospin-charge conjugate equation (IC equation), 

rW = (2/ 5 )(i|-m)$ (B.l) 

= (2/#)(i$ - m)^ ic , (B.2) 

from the left by a matrix of the general form \FrjF For each of these matrices, we 
will obtain three equations, the order of the list being: (a) subtract IC from non-IC 
equation, (b) add IC from non-IC equation, and (c) the result obtained by combining 
(a) and (b) to eliminate the bispinor with left-acting derivate operator, $ . Note that 
the equation we use for substitution will we written in terms of the aforementioned $ 
bispinor. 

Multiply by 1>: 

(g/2) J hv W hv = (i/2)(*0W - tf^tf ) - mJo (B.3) 
tf^tf = (B.4) 
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{g/2)J bv W bv = - mJ (B.5) 

Multiply by *r a : 

^rj> = ge bc J c v W bv - ^tM (B.6) 

(g/2) Jo v W w = (i/2)(*r a # - #V r tt) - mJ a (B.7) 

{g/2){JJ8 b + U c v e bc )W bv = MtJ* - mJ a (B.8) 
Multiply by 1>7 5 : 

*^7 5 tf = *7 5 #* + 2?miT (B.9) 

((?/2)i^W^ = (i/2)(* 7 5^* + ^75^) (B.10) 

{g/2)K bu W bv = itf 75 #tf - mfsT (B.ll) 
Multiply by "^r a75 : 

* ^75^* = "^i^O^a, - *75Tj* (B. 12) 

i^^X/W^ = (i/2)(* 75 rj* - ttV 75 T a tt) - ™K a (B.13) 

{g/2){K v 5 b + iX/e^)^ = i^rj* - m^ a (B.14) 
Multiply by 1>7 M : 

= -y^ftf - gS h i TW hu (B. 15) 

( ff /2) J 6 ^ = (i/2)(* 7 ^ - tt^tt) " mJ 0li (B.16) 

{g/2){j\ v - \S\ v )W bv = - mJ^ (B.17) 

Multiply by *r a 7 M : 

(<//2)( J„VV + 5 c /e a 6c )W 6l/ = (i/2)(*r a7/i # - ^V a7 ^) - mJ ¥ (B.18) 

*Vr„7„* = gWSeJ* - Sq/O W 6 „ - *r o7Ai ^ (B.19) 

(<7/2)(J <f/<f a 6 + iJ c 5/e a 6c - iS /5 a b + S c /e a bc )W bu = ittr^tt - mJ ¥ (B.20) 
Multiply by ^7 5 7m : 

( 5 /2)i^ 6 /^ = (i/2)(^ 757 ^ - * 7 57^) + mK 0li (B.21) 

^757m^ = -giK b W bfl - *757„^ (B.22) 

{g/2){K% v - i*S b ^)W bv = M-y^fa - mifo, (B.23) 
Multiply by *r 7 5 7„: 

*Vt„ 7 57m* = -i9(Ko&a% v + *S c /e a 6c )W^ - *r a75 7^* (B.24) 
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\{g/2){K c 5/e a bc - *S 0p »5 b )W bv = (i/2)(^r a757 ^ T albl /S>) - mK a , (B.25) 

(g/2)(K 5^5 a b + iK c 5^e/ c - i*S 0p »5 a b + *S c /e a bc )W bu = MT alblll $* - mK ap (B.26) 
Multiply by ^>a pe : 

= Va p Jy - (5/5/ - 5/5/)gJ b a W bv (B.27) 

(g/2)e p rK b a W bu = (i/2)(Va p JV + H^a^) - mS 0pe (B.28) 

(g/2)[e p rK b a + {{5/5/ - 5/5/)J b /W bu = Ma p J* - mS 0pe (B.29) 
Multiply by ~$T a a pe : 

(g/2)[e p rK J a b + i(5/5/ - 5/ 5/)ie a bc J c /W bv 

= (i/2)(Va p£ T a $y - ^a pe T/S>) - mS apt (B.30) 

^a pt r a ^> = g[{5/5/ - 5/5/)J J a b + e p /° e a bc K c /W bv - ~*o p( rM (B.31) 

(g/2)[e p r(K 0a 5 a b + iK CfJ e bc ) + i(5/5/ - 5/5/)(J 0a 5 a b + iJ CfJ e bc )}W bv 

= Ma^TafrS! - mS ape (B.32) 

Multiply by ^ 75 a pe : 

^ lb o p /S> = -U lb o pe $m - \e p / a gJ b c W bv (B.33) 

(g/2)i(5/5/ - 5/5/)K b a W bv = (i/2)(^ l5 a p J^ - ^ lf> a p ^>) - m*S 0pe (B.34) 

(g/2)[e p rJ b a + i(5/5/ - 5/5/)K b /W bu = i* 7 5<V#* - m*S 0pe (B.35) 
Multiply by ^r a75 cT pe : 

(g/2)[i(5/5/ - 5/5/)K 0a 5/ + \e p /° e a bc .J c /W bv 

= (i/2)(* 7 5<T pe r ^ - ^ lb a pe rj!) - m*S ape (B.36) 

^ 75 a pe r a * = g[(5/5/ - 5/ 5/)\e a bc K ca - ie p rJ J a b ]W bu - ^a pe rj^ (B.37) 

(g/2)[e p r(J 0a 5 a b + U ca e a bc ) + i(5/5/ - 5/5/)(K J a b + iK ca e/ c )]W bu 

= M l5 o- pe rJ^ - m*S ape (B.38) 



Appendix C. Fierz Identities for J-K Lorentz Vector Current Products 

Here follow the general Fierz identities for Lorentz vector currents J/, multiplied by 
the dual K/. There are four different types of Pauli term arrangements: a — b, — a, 
a — and — which indicate the % — 0, 1, 2, 3 index in Tj for the left and right terms in 
the product J/K/ . We treat the % — and i — a — b — 1,2,3 Pauli indices separately 
due to their slightly different algebraic properties. For each case, we write out the full 
Fierz expansion in spinor form, then expand Pauli and Dirac matrix products in terms 
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of irreducible elements of the Pauli and Dirac algebra using identities from appendix 
A. The final Fierz identities are written in the non-abelian current notation defined in 
Section 3. 

a — b case: 

J»K h v = ^7^(^)7,757^ 

= -(1/8) JiVw'YTaT'nV + (1/8) J la ^Yl a l v r a T l r b ^ 
-(1/16)^757^7^^^ 
+ (l/8)K ta ^Yl a l u r a r l r b ^ + (l/8)^WVr a rV^ 

= -(l/8)J„*75br " ^ u ¥ab + ie ab d r d ]V 

-(1/8) J c * l5 [rr - ia^][r a 5 b c + r b 5 a c - r c 5 ab - ie ab c }V 

+(1/8) Jo^t^V + ffV - <V - ^ x l5lx \[5 ab + ie ab d T d ]V 

+ (1/8) J^^rT-f + V au Y - rf v l° ~ ^ x l5lx }[rJ b c + r b 5 a c 

-r c 5 ab - ie ab c ]ty 
-( 1/16) S 0ae ^ 5 [i^rf" - \rf v rf e + ifo*" - rf v - e aeu ^ 5 

+ie^ x l5 ^ x ][5 ab + ie ab d r d ]^ 

-(l/16)5 CCTe ^7 5 [i77 e V CT - i?f V e + rFc 1 " ~ ~ ^""75 

+ie^ x l5 a\][r a 5 b c + r b 5 a c - r c 5 ab - ie ab c ]V 

+ (l/8)K 0a ^[v f " a l u + V au l" ~ V^l a ~ ie^ A 7 5 7A][<U + ieJu]* 
+(l/S)K c Mv tMT r + yf v r - v" u l a - ie^ x l5lx ][r a 5 b c + r b 5 a c 

-r c 5 ab - ie ab c ]^ 
+(l/8)K *[rr - ia^][5 ab + ie ab d r d }V 
+(l/8)K c *[rr - ia^][r a 5 b c + r b 5 a c - r c 5 ab - ie ab c ]V 

J/K b » = (1/4) [U a *Sr + iJ b *S a ^ - \K a Sr - iK b S/ v + J*K h v 

+J a »K b » + J b "K a v + J b v K a » - J^Kf-fT - JtoKSrT 

+5 ab (iJ Q *Sr - iJ c *S c ^ - \K Q S<T + iK c S^ + J Q »K " 

+J v K > i - J/K c » - J»K C » - JtoKfrT + Jc*K°*rr)] 

+ (l/4)e ab c [-U K c r ] ^ + iK J cV ^ + J ca J ox e^ x + K ca K 0X e^ x 

+(l/2)i(-S ^*Sr - S Q \*Sr + S C W + S C \*S^)\. (C.l) 

a — case: 

JSKo" = ^7^(^)757^ 

= -(1/8)^757^7^^^ + (1/8) JirV'ft'f'f'fTaT'* 

-(1/16)5^757^^7^^^ 
+ (l/8)K^7 M 7 CT 7^ a r^ + (l/8)K^Yr a r 1 ^ 
= -(1/8) J *75br - i^]r a * - (1/8) J c * 75 [7r - ^"IP-" + ie^]* 
+ (1/8)^75^^7" + ?T7 M - <V - ie^ A 7 5 7AK^ 
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+(1/8) J^TsbTT" + V au Y ~ vTl - ie^ A 757A]K c + ie a cd r d }^ 
-(l/16)5o CTe *7 5 [i77 ei/ 77^ - \rTfT + ffa^ - rf v a^ - 6^75 
+ie^ x l5 a\]r a ^ 

-(l/16)5 CCTe *7 5 [i77 e V CT - < V + rf v o^ - rf v a^ - 

+ie^ x l5 o-\][5 a c + ie a cd T d ]V 
+ (l/8)K 0a ^[r]^Y + rf v r - rT-f - ^ x ^ lx }T a ^ 
+(l/8)K ca V[r)^Y + - - ie^ A 757A][5a C + ie a cd r d ]V 

+ (l/8)K ^[r]^ - icr^]T a * 
+(l/8)K c *[rr - io-^][5 a c + ie a cd r d ]V 

J*Kf = (1/4) [iJJSr + U a *Sr - iK Sr - iK a S<T + 

+Jo u K a » + J a »K u + J a v K^ - JtoKSrT - J„ CT X V1 
-(l/A)e a cd [U c K dV ^ + (l/2)J ca J dX e^ x + (l/2)K ca K dX e^ x 
-(l/2)i(S/ a *Sr + S d \*S°»)\ (C.2) 

— a case: 

J ^/ = W(^K757^ 

= -(1/8) J^7 5 7^TV a * + (1/8) J^s-f-fYr^Jf 
-(l/16)^ e ^757^ CTe 7 I/ rV a * 
+ (l/8)^^ 7 CT 7^V a ^ + (1/8)^7^7"^* 
= -(1/8) Jo* 75 [»r " i^]r a * - (l/8)J c tt 75 br - i<7"1P« c - ie/^]* 
+ (1/8)^75^^7" + " <V - ie^ A 7 57A]r a ^ 

+ (1/8)7^75^^7" + V av Y - vT-f - ie^ x l5lx \[5 a c - ie a cd T d ]V 

-(l/lQ)S ^l5[iv eu V^ ~ iv au V^ + V eu ^ a ~ V au ^ e ~ e CT£ ^75 
+ie^ x lb a\]r a m 

-(l/16)5 CCTe *7 5 [i77 e V CT - 'rf v rf e + V eu ^ a ~ V™^ e - ^^75 
+ie™ x l^ x ][5 a c -it a cd r d }m 

+ (1/8)7^(^7" + TT-f ~ ~ ie^ A 7 57A]r a ^ 

+ (l/8)K ca ^[ V ^Y + V™Y - <V - ie^ A 757A][5a C - ie a cd r d ]^ 

+(l/8)K *[rr ~ 

+(l/8)K c *[rr - ia^][5 a c - ie a cd r d ]^ 

J »K a » = (1/4) [i J *Sr + U a *S^ - \K Sr - iK a Sr + JfK a v 

+J v K a » + J^Ko" + J/tfo" - 3^K a a rT ~ J^rT] 
+ (l/4)e a cd [U c K d r ] ^ + (1/2) J ca J dX e^ x + (1/2)7^7^^ 
-(1/2)1(5//^ + S d \*S°»)\ (C.3) 

— case: 

JfK Q v = ^7^(^)757^ 
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= -(1/8) J^fcYVV* + (1/8) J iCT * 7 57 / VYV* 
-(1/16)^757^7^ 
+ (l/8)^ (T '*7 At 7 CT 7^ i * + (1/8)W7^ ! * 

= -(1/8) Jo*T5[^ - - (1/8)J C *7 5 [77^ - ur^r"* 

+ (1/8) Jo^iv^r + V au r ~ rT-f ~ ie^ A 7s7A]^ 
+ (1/8)7^75^^7" + rT-f ~ Tf v l° ~ ie^ A 7 57A]r c ^ 
-{l/16)S 0m V'y s [iri a 'ri tMT - + r] ev a^ - if a** - e aeu ^ 5 

-(l/lQ)S CCTe ^ 5 [ir] eu r]^ - irTrf* + rf v - rfa^ - e atu ^ 5 

+(1/8)K ^[ V ^Y + rfV - rT'f ~ ie^S™]* 
+ (l/8)K c JJ[7 ] fieT Y + f / - r/'V - ie^ A 7 5 7 A ]r c ^ 
+(l/8) J RT ^[< iy - 10^"]* + (l/8)i^ c ¥[^ - icr^]r c * 

J^K Q V = (1/4) [iJo*^^ + U*S cta/ - iKoSo^ - iK c S Cflu + J^K V + Jo^o" 

+J/K CU + J/iC* - J^KJrf v - J„K w rT] (C4) 

Appendix D. Fierz Identities for Rank-2 Skew Tensor Currents 

This appendix contains a more detailed version of the derivation of the Fierz identities 
for So'"', *S M ", S a ^ u and *S a M ". Let us first derive expressions for S ^ and *S M ", using 
the Fierz identities for JK vector current products derived in appendix C. Consider the 
following expression, obtained from the Fierz identity (jC.4p and subtracting the Pauli 
trace of (IC.lj) : 

WKo" - J a u K a » = U Q *Sr - iK Sr - (1/2) J »K V - (1/2) Jq-W + {l/2)J a »K w 

+(1/2) J a v K«» + (1/2) JtoKfrT ~ (1/2) J^V"- (D.l) 

Note that we have used the Fierz forms that have no stand-alone Pauli vector triplet 
indices a = 1,2,3. Similarly, using the same Fierz identities we can also form 

JaW - J V K » = U *So" v - iK S " v + (1/2) J »K " + (1/2) J ^o" - (1/2) J a »K av 

-(1/2) J/iT" - (1/2) J^Ko'rT + (1/2) J^V"- (D.2) 

Adding these two equations, we get 

j»K iv - J?K* = 2i(J *S(T - K S ^) (D.3) 

which is antisymmetric in /i, v. Another Fierz expression which is antisymmetric in /i, u, 
is 

e^ vpk Jq p Kq k = (l/A)e^ pR [U *S 0pK + iJ a *S a pK — iK So P K — iK a S a pK + Jq p Kq k + Jo K K 0p 
+J ap K\ + J aK K% - J aK a Vp K - Ja„K a °ii pK \ 
= -(1/8) (e^e^) J S 5r - (l/8)(e^% K5T ) J a S a5r - (1/2)K *S ^ 
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-(l/2)K a *S a ^ 

= (1/2) Jo + (1/2) J a S a » v - {l/2)K Q *S<r - (l/2)K a *S aflu (D.4) 

where we have first canceled out all of the terms from flC4l) symmetric in /i, u, then 
applied the Levi-Civita double contraction identity 

6^%^ = 2((JV5V - <*V*V) (D.5) 

in addition to the definition (|26|) to convert the S terms to *S, and vice-versa. The 
other term we can contract with the Levi-Civita tensor, with no stand-alone Pauli 
vector triplet indices is the trace over the Pauli index a. So after (again) canceling all 
symmetric terms from (IC.lj) and applying the Levi-Civita contraction identity, we get 

e» vpK J ap K\ = (3/2)J S<r - (1/2) J a S a ^ - (3/2)JW + (l/2)K a *S a ^ (D.6) 
Taking the sum of (1D.4j) and (1D.6j) . we obtain 

e^J ip K\ = 2{J Q Sr ~ KfStT). (D.7) 
Now take the combination of (1D.3I) and (ID. 71) 
J e^ K J ip K\ - iK Q (JSK iv - J?K^) 

= 2{JlSr - J Q K *Sr + J Q K *Sr - KlSr). (D.8) 

Canceling the middle two terms on the right-hand side and rearranging, we end up with 
one of our identities 

Sr = (1/2)(J 2 - KlY l [J Q e^J ip K\ - iKoiJfK" - JfK*")]. (D.9) 

By comparison, the abelian version of the identity is |llj . |15] , 

sT = {a 2 - cu 2 )- l [ae^ Jp k K - iu{fk v - fk")], (D.10) 

where we have used an alternate definition for the pseudoscalar, u = V^/s^S as opposed 
to the often used i/ji^ifj. Now to calculate the dual, take an alternate combination of 
( 1D~3|) and ( lD~7j) 

K e^J ip K\ - U Q (JfK iv - J?IC») 

= 2{J K Q Sr - K 2 *S ^ + J 2 Q *Sr - J K S n (D.H) 
which immediately leads to 

"ScT = (1/2)(J 2 - KlY\K^ vpK J ip K\ - U Q {JfK iv - J»K^)\ (D.12) 

Let us now derive an expression for S a ^ v . We start by forming an expression from the 
JK Lorentz vector current product Fierz identities, with a single Pauli vector triplet 
index, (E^]) and f lU3|) . 

j/jv + j^K a v = {i/2)[xj Q *sr + yJa*sr - \K sr - \x a s<r + wks 

+J u K a » + J a »K u + J a u hV - JoaKa'rT - JarKfrT], (D.13) 

which eliminates the Pauli Levi-Civita contracted terms that are in the separate a — 
and — a cases, as they are identical and of opposite sign in each. Subtracting (ID. 131) 
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from itself, but with the /z, v terms flipped, will cancel all of the terms symmetric in 
[M, v. 

(J a »K » + J "K a u ) - {J a v K Q » + J u K a ») 

= ujsr + u a *sr - iK Q sr - (d.u) 

Another way to get rid of the terms symmetric in fi, v is to contract (ID. 131) with e^ vpK : 
e^ K (J ap K 0K + J 0p K aK ) = (l/2)e^ K [U *S apK + U a *S 0pK - iK S apK - iK a S 0pK ] 
= -(l/A)(e^e pK5r )J S a ST - (l/A)(e^% KST )J a S Sr - K *S a ^ 

= -{1/2){8» T 8 V 5 - 5»s5\)JoS a ST - (1/2)(5» T 5 V S - 5^5\)J a S Q 6r 

-K *S a ^-K a *S » u , (D.15) 
where we have used the Levi-Civita tensor double contraction ( ID. 51) . Simplifying gives 
e^ K (J ap K 0K + J 0p K aK ) = W + J a Sr - K Q *Sr - K a *Sr- (D.16) 
Now take the combination of ( ID. 141) and ( ID. 161) : 

J e^ K (J ap K 0K + J 0p K aK ) - iK [{J/K^ + JfK a v ) - (J a v K^ + J U K^)} 

— J^Sa^ + Jo JaSo^ — JqKq*SJ 1U — JqK^Sq^ + Kq Jq*S ci ,J ' V + KqJ^Sq^ 

-Kisr - K Q K a Sr 

= (J 2 - Kl)Sr + Wa - K Q K a )Sr + (K J a - J K a )*S<r. (D.17) 

Let us rewrite the second and third terms on the right-hand side, by substituting the 
So^ and *S t ™ identities, flD~9l) and flD~T2l) : 

(J J a - K K a )S ^ = (1/2)(J 2 - Kl)~ l [{JlJ a - J K K a y upK J ip K\ 

+i(K 2 K a - J K J a )(J^K iv - JfK*)]. (D.18) 
Likewise, for the dual we have: 

(K J a - J Q K a fSr = (1/2)(J 2 - Kly^Kpa - J K K a )e^ K J ip K\ 

+i(J 2 K a - J Q K J a ){JfK w - JfK**)]. (D.19) 
Summing these two equations together gives 

J 2 n + Kl 



(J J a - K G K a )Sr + (K J a - J K a )*S 



2(J 2 - K 2 ) 
J K 



T 2 - K 2 



[J a e^J ip K\ + iK a {JfK w - J»K^)} 
[K a e^ pK J ip K\ + U a (JfK lv - JfK^)]. (D.20) 



Substituting this into ( 1D.17I) and rearranging, we finally get 

S a ^ = (J 2 - Kl)- l {J e^(J ap K 0li + J 0p K aK ) 

-iK [(JSK " + WK a ») - (J a »K » + J »K a »)}} 
J2 + K 2 



'2(J 2 -K 2 ) 2 



[J a e pupK J ip K\ + \K a {JfK™ - JSK i n] 



+ (T2 J ° K L, 2 [Ka^ pK J ip K\ + iUJfK™ - JfK**)]. (D.21) 
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Now let us derive the dual, *S a ^ u by considering the alternate combination: 
K e^ K (J ap K 0K + J 0p K aK ) - iJ Q [WK Q v + J »K a v ) - {J a v K^ + J Q v K/)\ 

= j K Q sr + K Q j a sr - Ki*sr - K G K a *sr + 4*sr + wsr 
= (j 2 - Kiys/» + (K Q j a - j Q K a )sr + (j j a - K K a ysr. (D.22) 

Again, we rewrite the second and third terms, by substituting identities (1D.9|) and 
(1IXI2]) : 

(K J a - J Q K a )S<T = (1/2)(J 2 - Kl)-%J K J a - JlK a )e^J ip K\ 

+i(J K K a - KlJ a )(JfK iv - JfK^)]. (D.23) 

The dual term is 

(J J a - K Q K a fSr = (1/2)(J 2 - Kl)-%J,K J a - K 2 K a )e^ pK J ip K\ 

+i(J K K a - JlJ a ){JfK w - JfK**)]. (D.24) 

Summing the two: 

oJa — K K a )*So fJ ' 1 ' 

J K 



(K J a - J Q K a )Sr + (J J a - K K a yS 



J o ^0 



' J '' A '' ) \K a ^ p \J ip K\ + iJ a {JfI& - JfK*)]. (D.25) 



2( J 2 - 

Finally, after substituting this into (ID.22[) and rearranging, we get: 
*S a ^ = (J 2 - KD^iK^^JapK^ + J 0p K aK ) 



J 2 o + K% 



J " /V " :[J a e^ K J ip iT K + iK a {JfK iv - J t v K^)\. (D.26) 



Comparing this identity with the identity for S a ^ v ', flD.21j) . we can make the observations 



that the first terms are the same in each, but with the Jo and Kq terms interchanged 
(similarly to the Sq' 11 ' and *So fll/ comparison). Likewise, the second and third terms of 
both identities are the same, but with the J a and K a terms interchanged, along with the 
signs, which are flipped with respect to each other. 
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